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Abstract 

In this paper, we study phase transitions in asymmetrical fermion superfluids. In this scenario, 
the candidates to form pair are particles with mismatched masses and chemical potentials. We 
derive an expression for the critical temperature in terms of the gap and masses (or chemical 
potentials) when the constraint of equal Fermi surfaces m a ^ a = mft/ij is imposed. 
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The advent of new techniques to deal with ultracold fermionic atoms has motivated the 
interest on its theoretical investigation. One of the most intriguing systems that have been 
investigated in last few years are the asymmetric ones. In these asymmetrical Fermi systems 
the masses and densities or chemical potentials of the two species that will possibly form pair 
are unequal. Candidates for the ground states of these cold, dilute quantum systems have 
been proposed: the Sarma phase [jj, the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) phase 0, 
the breached-pair superfluidity (BP) phase and the mixed phase (MP) 0,0|- Pairing 

in asymmetrical fermionic systems are predict to have applications from the explanation of 
chirality among amino acids [H, 0] to laser atomic traps experiments [T3| • An immediate and 
important question that naturally arises in this fascinating field is: How does the critical 
temperature of an asymmetric system depend on its asymmetry? 

The critical temperature is the temperature at which the gap vanishes signing that the 
system is in the normal state with unpaired particles. If the mass and density asymmetries 
are large, the attraction is weak and we expect that in this asymmetrical conditions the 
critical temperature is smaller than that of the symmetric system. If the asymmetries 
increase more, even a small amount of heat is enough to break the Cooper pairs and disrupt 
superfluidity. 

In this paper we investigate two species systems, and answer this question for two possi- 
bilities for the particle's masses and fixed chemical potentials asymmetries. We show that 
in a fermionic system constrained to m a fi a = mbUb, with m a ^ mb and \x a ^ fib, the critical 
temperature is a slightly decreasing function of the asymmetry. We find a generalization for 
the expression relating the critical temperature and the zero temperature gap parameter. 



1. The Model 

As mentioned earlier, several candidates forpairing of unequal Fermi surfaces have been 
proposed in the literature. The Sarma phase 0] and the prior version of the BP phase 
are unstable (ESQ; an d w iU n °t be considered. Regarding the LOFF state, since this 
phase can exist only in a very narrow window of asymmetry for the Fermi surfaces, we do 
not consider this realization of superfluidity in the present calculations. The investigation of 
the critical temperature in the MP (with fixed particle densities) 0,0]' anc ^ m the clarified 
version of the BP phase 1 (with fixed chemical potentials) 0, must be considered and will 
be discussed elsewhere. 

We employ the Bardeen-Cooper-Schrieffer (BCS) model to derive the Ginzburg-Landau 
(GL) theory, following the modern formulation developed by Sakita [l]]]. For the sake of 
completeness, we also derive the gap equation for an asymmetrical fermionic system at finite 
temperature by the variational method. The derivation is shown in the appendix, and the 
results agree, as it should. 

We consider an asymmetrical nonrelativistic dilute system of a and b fermion species 2 , 
having masses m a and m&, with Fermi momentum P 3 F = ^Jlm^fi^ j = a, b. 



1 This version takes into account the momentum structure of the interaction and large mass asymmetry. 

2 In (relativistic) quark matter, the stud y o f the critical temperature from the point of view of the GL 
approach, has been carried out in Refs. Ill El. 
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Let us begin with the partition function 

J D[^\D[tfjel- s W\ (1) 
where S(ip, tft) = — J dr J dxL, with L being the BCS Lagrangian 

L =J2 ^ [~ d - + ^ + A ^ + 9^W b {x)^a{x)^{x), (2) 

i=a,b 

where g > 0. Introducing auxiliary fields via the Hubbard- Stratonovich transformation, we 
find an effective action, in which the mean field BCS Lagrangian is expressed by 

L MF = ^il>l(x)(-d T + ^ + lHjil>i(x)+ (3) 

i=a,b ^ 1 ' 

Afx)| 2 



A(x)ipl(x)i/jl(x) + A*(x)i/j b (x)ipa(x) + 

y 

Introducing a source for the A(x) field, we write 

Z = N J D[lf)]D[^]D[A]D[A^] e JdrJdx[L MF+ rA( X ) +j A*(x)] ^ 

The partition function is Z = Z\j,j*]j=j*=o and the generating functional for the connected 
Green's functions is defined as 

W\jJ*]=hxZ\jJ% (5) 
The Legendre transformation of W\j, j*] is 

r[A, A*] = I d A x(fA(x) +jA*(x))-W[j,f}. (6) 

From all one-loop diagrams contributing to T[A, A*], we evaluate only the one which gives 
a contribution for the A 2 term. It is shown in Fig. ((TJ. 

2. The one loop correction to the BCS gap parameter 
The effective action up to one-loop is 

T[A, A*] = a\ A| 2 + (3\ A| 4 - C A*^V 2 A, (7) 

where M = mamb is the reduced mass, a — - — A, and A is the momentum independent 

m a +m b g ' L 

contribution from S(g = 0, q), which is given by 

sfa = o, a = g (faw _ ,. W)( U _ ^ m = ( 8) 




FIG. 1: The one-loop diagram contribution to T[A,A*] with two external lines. 



(|7j) is valid only in the vicinity of a secondj^rder phase 



where u n = (2n + 1)tiT, oj a {k) - 
It is worth noting that Eq. 
transition, since it is based on the assumption that the gap parameter is small [14j. While it 



is important that A be small so that higher order terms can be neglected, it is also important 
that the fermions have a finite gap so that they can be properly integrated out to obtain the 
effective potential. It is this latter condition that renders the present analysis appropriate 
only when the Fermi surfaces are not mismatched. 

After the frequency summation, employing the imaginary time formalism of finite tem- 
perature field theory, we obtain 



.4 



d 3 k 



(2vr) 3 uj + u' 
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(9) 



where, for short, we have defined to = cu a (k), uo' = u b (k) = ^ — and j3 = where k B 
is Boltzmann's constant that we will set equal to one. Changing the variable of integration 
from k to e, we get 



.4 



Wc de 
e 

^ de 



1 



r c de 
P(0) / - 
Jo £ 



e f3u>{e) _j_ I 
1 



e Pu><(e) 
1 



(10) 



e Puj(e) _|_ I e Pu>'{e) _|_ J 



where p(0) = is defined as the density of states at the Fermi level, with &p = \J2M^i 
being the "average" Fermi surface and \i = \i a + We have also defined W C (A) = 
a; (A) + a;' (A), where A is the cutoff in the momentum integral. The energies are given by 



M M 

u(e) = — e + n n a , 

m a m a 

M M 
uj (e) = — e + /I fi b . 



(11) 
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After simple algebra, Eq. ([lUp can be written as 



.4 



p(0) f Wc de 



tanh(^ + tanl/^ 



If we define 



M ( M 

V = H ~ Ha = - A 4 - 

m a \m b 



then we can write Eq. (|12|) with the aid of Eqs. and (jlHj) as 



(12) 
(13) 



,4 



2 \y e V2m a 2 1 
P(0) 



Wc cfe . f (3 M (3 
tanh ( — — e 7] 



tanh (x + a) 



tanh (x — a] 



x 



2 rrib 



(14) 



where we have defined 

(3M (3M (3 

\ a = - — W c , A fe = - — W C and a = -r]. 
2 m a 2 nib 2 

If we take > m a , then Xf, < X a and Eq. $T^i can be written as 



(15) 



.4 



P(0) 



Xa dx , , . /" Ao 

— tanh [x + a) + — tanh (x — a) 

o x io 



Aa dx 

— tanh ( x — a 

b x 



(16) 



We can solve the first two integrals in the r.h.s. of Eq. (|16|) employing the residue theorem 
and the last one can be easily solved if we observe that A a and \ correspond to the regions 
where tanh(x + a) ~ tanh(x — a) ~ 1. Thus, 



(17) 



where ^(a) = \1>(| + — ) + \1>(| — — ) with \1/ being the digamma function, defined as $(z) = 
wjy, where z is a complex number with a positive real component, T is the gamma function, 
and r" is the derivative of the gamma function. We also have that ^(0) = —27 — 41n(2), 
where 7 is the Euler's constant. 

The critical temperature is the solution of the equation 



1 

a = - - A = 0. 
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(18) 



Then we write 



In (j3a- 



u D \ 1 
gp(0) *~ 7i ) = 2 



S*(a) 



(19) 



where a 



M 



y/m a m h 

used the fact that = 2uj£> 



is a dimensionless parameter reflecting the mass asymmetry, and we have 
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The BCS gap in the weak coupling limit, p(0)g « 1, is given by A = 2uj£) e~ 1//p (°) 9 . 
Since the gap which minimizes the free-energy of the asymmetric system retains the same 
size as in the symmetric case (An) until a value for the asymmetry where the pairing is not 
afforded any more |E 0, EE IE 13] > we rewrite Eq. (fT§j) as 



T c = G -pLe-^\ (20) 

where a c = y?7 = ^ ^m^^" ■ We evaluate Eq. ()20|) in the two possible configurations 
for the particles masses and chemical potentials constrained to P F = P F , which are the 
situations encountered when the fermions are fully gapped. In these cases we can obtain 
analytical solutions for the critical temperature, due to the simple form of the term jF(a c = 0), 
as showed below. 

In the investigation of the phase transition when the Fermi surfaces are mismatched, one 
needs to compare the thermodynamic potentials of the superffuid and normal states. In 



fact, as showed in p, [D, [15|, |12| by the behavior of the free-energy as a function of the gap 
for several asymmetries in the chemical potentials, an asymmetrical fermion system stays 
in the superfluid phase until a maximum value for the difference in chemical potentials is 
reached. After this maximum value, there is a first order phase transition to the normal 
phase. 



3. Equal Chemical Potentials and Masses 

This case configures the symmetric system, whose critical temperature is recovered for 
m a = rrif,, Ha = /J-b, resulting 3(a c = 0) = — 21n(4e 7 ), and a = |, giving the well known 
BCS result 

T c = — A = T° ym . (21) 

7T 

We notice that both the zero temperature gap parameter An and the symmetric critical 
temperature T^ ym depend on the product p{0)g, but not their ratio. 



4- Equal Fermi Surfaces with Mismatched Chemical Potentials and Masses 

This situation is also characterized by P F = P F , which implies m a p a = m^Pbi however 
with m a 7^ m b , p a ^ p b . This is achieved by setting a c = in Eq. (|20|). yielding 

T C (P F = P b F ) = 2a— A = 2aT s c ym . (22) 

7T 

We note that if we set m a = nib in the equation above (which implies a — |), then we 
have the symmetric Fermi gas, since we would also have p a = pb, due to the constraint 
f^aPa — rribPb- We can observe from Eq. (|22|) that the critical temperature for the system 
constrained to m a p a = rribPb goes with 2^j^- T^ ym for nib greater than m a and approaches 

zero for nib » m a- This shows that the pair formation is disfavored for very large mass 
asymmetry. The same conclusion has been found for the case of fixed number of particles 0] . 
In Fig. (J2J we show the ratio T c /A as a function of nib/m a . As one can see, T c /A is a 
smooth function of the mass asymmetry, and goes to zero for nib » ma- 
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FIG. 2: T c /Ao as a function ofmb/m a for a system constrained to m a (i a = mbfij,. 

A mass ratio in the order of, or larger than, m^/rria = 50, as used in Ref. |1| would, 
obviously, have a smaller critical temperature than that of the symmetric gas. Although 
we did not consider in our calculations the momentum structure of the interaction em- 
ployed in Ref. [Hf, we believe that our results would persist, at least qualitatively, after this 
consideration. 



5. Conclusions 



We have calculated the critical temperature of an asymmetrical Fermi system in two 
configurations for the masses and chemical potentials of the two species that form Cooper 
pairs. Among those cases we have investigated, the one constrained to vn a [i a = rribfib, with 
wia 7^ Tt6 ; Ha 7^ ^b-, is particularly interesting, for which we found a generalization for the 
expression relating the critical temperature and the zero temperature gap Ao- Namely: 

A 



and constitutes an universal constant, for 



. jlI. w 1 76^ where a = = VJL ^ i 

T c e~> 2a 2cr ' m a +m b [ta+Hb 

given m a and m& (or fi a and independent of g and p(0). Another remarkable feature 
of this result is its independence of any cutoff parameter. This is because Eq. PT)|) is quite 
insensitive to the regularization procedure. 

We believe that the results achieved in this work could, in principle, be tested experi- 
mentall y in , for example, experiments involving G Li or m K in atomic traps 0, H EI 



231 . |2J, |25|, |2a, |27j, |2jj |29|, |30j, |3JJ, |32j. Evidences of superfluidity in these s yste ms were 




observed both microscopically, observing the pairing of fermionic atoms 12^ 
macroscopically, due to anisotropic expansions [2j|, collective excitations [2^ 
heat capacity measurements jjyj. In these experiments, the strength of the pairing inter- 
action can be controlled by an applied magnetic field, for instance. Also, the density, the 
number of each species and the trapping potential can be altered. Thus, the species on the 
trap differ by theirs spin/pseudo-spin projections, and also by theirs densities, characterizing 
the asymmetry of the system. 
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Current experiments 0, 3^, 31, 33] produce temperatures down to about 0.057V, where 
Tp is the Fermi temperature for a noninteracting gas with the same number of atoms and 
trap conditions as the experiment, typically of order of fiK. However, a weakly interact- 
ing Fermi gas requires much lower T to achieve superfluidity. For the conditions of these 
experiments, the mean field approximation with an interaction energy proportional to the 
scattering length is not valid. However, the mean field approximation with a unitary limit 
appears approximately valid, furnishing a good agreement with predictions of the collective 
frequencies, and a very good agreement on the transition temperature (33[. Thus, even 
when the measurements are done in strongly interacting Fermi gases, mean field theory had 
qualitatively explained the behavior of these systems, and we expect that our weak coupling 
mean field BCS results should also be valid qualitatively. In particular, we believe that the 
numerical factor 1.76^- can be tested experimentally, provided the mass asymmetry is not 
so large. 

When the Fermi surfaces are mismatched the phase transitions will be of first order and 
the present formalism is insufficient to determine the critical temperature. We will come 
back to this issue soon [34j|. For these cases m a \i a ^ m b fi b , and ^ will not be an universal 
number. This "non universality" is also manifested in the BP at finite temperature [35| 
and in (dense) neutral quark matter jH, E^E^, which is essentially asymmetrical 3 . Still 
in quark matter, in Ref. [39] was developed a systematic method of QCD expansion of 
the transition temperature, motivated by the non-BCS scaling of the gap parameter with 
coupling. In this work, the relation between the zero temperature energy gap and the critical 
temperature has a non BCS form too. 

We find that large mass and chemical potentials asymmetries lower the critical tempera- 
tures substantially and compromise the stability of the system. This happens because any 
small thermal excitation breaks the (weakly bound due to the large asymmetry) pairs, and 
destroys superfluidity. 
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6. Appendix: Solution by the Variational Method 

We now derive the gap equation at finite temperature for an asymmetrical fermion sys- 
tem, in order to determine the critical temperature. We follow the usual derivation of the 
textbooks, however extending the analysis for the asymmetrical systems we are investigat- 
ing. Let us define fk as the probability of an a particle with momentum k is excited, and 
similarly as the probability of a b particle with momentum — k is excited. Then, the 



3 The up, down and strange quark Fermi surfaces, which are candidates to form color superconductivity are 
mismatched. 



s 



entropy for an asymmetrical fermion gas is found to be 



S = ~ E U' k ln ^ fc ) + ( X " /*) ln ( X - /*) + 9k H9k) + (1 - 9k) ln(l • (23) 

fe 

The free energy or thermodynamic potential is written as 

F = E — TS, (24) 

where E is the internal energy 

E = E WK 1 ~ & ~ ^ + /*! + * ~ & ~ ^ + 9k]} (25) 

-gy^^fc'Vfc'^fcVfcti - Sh - fl , fc)(i - fk> - 9k')- 

k,k> 

Here we have defined the particles energies relative to their Fermi surfaces in terms of our 
previous definitions e£ = to = ^ an d e t = CJ ' = 2^ ~~ A* 6- From the minimizations 

f = °- (26 » 

ogk 
5F 

0, 



we find, respectively, 



A = 1/(^ + 1), (27) 
^ = l/(e^ + l), (28) 



^=K i+ tj' (29) 

where £^' b = ±e^ + are the quasiparticle excitations, with El = e^ 2 + A 2 (T) and 
e\ = 6fcI ^ €fc . In the definition of £^' 6 we have also defined 

A(T) = gJ2UkV k (l-fk-g k ). (30) 

k 

Since V k 2 = 1 — U k , then = 7^-, and the equation above can be written as 

l= 9Y,w {1 ~ h ~ 9h) - (31) 

k 
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At the critical temperature A = in Eq. (|31|). and we have 



k ~ k 



which is equation (|18|). The solution for T c from (|32j) follow as in the body of the paper. 
Important remarks are now in order: 

1. Although Eqs. (pTTj) and (}3"2j) require regularization, the regulator dependence cancels 
from the result P0"j) . 

2. We have obtained an equivalence between two approximations to identify the zero 
temperature quasiparticle excitations. Introducing the temperature via the variational 
method, minimizing the thermodynamic potential with respect to the excitations proba- 
bilities in thermal equilibrium and then taking the zero temperature limit is equivalent to 
diagonalize the (zero temperature) mean field Hamiltonian (as done in Refs. 0, 0|) and 
obtain the excitations energies of the fermion quasiparticles. 

3. If the fermions are fully gapped, as is the case here, the gap parameter depends on the 
asymmetries only at finite temperature (0 < T < T c ). At zero temperature the excitations 
probabilities vanish, and the gap depends only on "averages" , through et. When there are 
gapless excitations, the gap depends on the asymmetries even at T = 0,0, IE Q- 
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